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ABSTRACT 

We report Hersche^\ SPIRE (250, 350, and 500 /im) detections of 32 quasars with redshifts 0.5 < z < 3.6 from 
the Herschel Multi-tiered Extragalactic Survey (HerMES0)- These sources are from a MIPS 24 /im flux-limited 
sample of 326 quasars in the Lockman Hole Field. The extensive multi-wavelength data available in the field 
permit construction of the rest-frame Spectral Energy Distributions (SEDs) from ultraviolet to the mid-infrared 
for all sources, and to the far-infrared (FIR) for the 32 objects. Most quasars with Herschel FIR detections show 
dust temperatures in the range of 25 K to 60 K, with a mean of 34K. The FIR luminosities range from 10 11 3 
to 10 13 5 L Q , qualifying most of their hosts as ultra- or hyper-luminous infrared galaxies. These FIR-detected 
quasars may represent a dust-rich population, but with lower redshifts and fainter luminosities than quasars 
observed at ~ 1 mm. However, their FIR properties cannot be predicted from shorter wavelengths (0.3-20 /im, 
rest-frame), and the bolometric luminosities derived using the 5 100 A index may be underestimated for these 
FIR-detected quasars. Regardless of redshift, we observed a decline in the relative strength of FIR luminosities 
for quasars with higher near-infrared luminosities. 

Subject headings: galaxies: active, nuclei, starburst - infrared: galaxies - quasars: general 



1. INTRODUCTION 

The far-infrared (FIR) properties of quasars are critical for 
our understanding of active galaxy evolution, as they probe 
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the possible evolutionary connection between star forma- 
tion and black h ole (BH) accretion (e.g. ISanders et al.lll988t 
iLutz et all [2007). Both starbursts and active galactic nuclei 
(AGN) could contribute to the rest-frame FIR emissions at 
various redshif ts, although starbursts a re considered to b e 
dominant (e.g. Rowan-Robin sonlll995t Pfrichas et aLll2009h . 
Other models have successfully explained FIR luminosity as 
originating from direct AGN heating, where cool dust at large 
distances from the AGN reside in a warped disk/torus (e.g. 
ISanders et aflll989l; iHaas et al.ll2003l) . Spectral information 
can be used to break this degeneracy. Polycyclic Aromatic 
Hydrocarbon (PAH) features, for instance, are us ed to indi- 
cate star formation activity (e.g. ILutz et al.1 120071) . Spitzer 
spectroscopic studies detected PAH emission in some of the 
Palomar-Green (PG) and Sloan Digital Sky Survey (SDSS) 
quasars, where strong star formation coexists with quasars, 
and is re sponsible for an average of 60% of the FIR emis- 
sion (e.g.lHao et alJ2005tlNetzer et al .l2007l:lLutz etaLll2008h 
IVeilleux et al.ll200StlShi et alJl2009t) . 

Rest-frame FIR SED studies provide powerful con- 
straints on the star formation in quasars. Despite 
large dispersions at different wavelengths for optically 
bright, unobscured quasars, their mean SEDs show sur- 
prising unifo rmity over redshift, luminosity, and Edding- 
ton ratio (e. g. Elvis et al. 1994, hereafte r E94l [Richards et alj 
12006. herea fter R06; Hao et al.ll201 ll) . However, the SEDs at 
rest-frame A > 40 /im were poorly defined at high redshifts 
(IRAS 100 /im was the longest wavelength for E94 and R06). 
Several groups have tri ed to address this FIR gap with various 
samp le selections (e.g. lPapovich et ai1 l2006: Kartalte pe et alJ 
2010); others have taken advantage of (sub)millimeter ob- 
servations, e.g., with the Infrared Space Observatory (ISO), 
the Institut de Radioastronomie Millimetrique (IRAM) tele- 
scope, and the Sub millimetre Common-User Bolometer Ar- 
ray (SCUBA) (e.g.lOmont et al.l200lLl200llHaas et alEool 
Priddey et al. 2003a). However, a single photometric point in 
the rest-frame FIR, as is the case for most (sub)mm quasar 
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studies, does not strongly constrain the dust temperature dis- 
tribution for these quasars. The sample sizes for (sub)mm 
quasars with enough photometric points that allow FIR SED 
studies, on the other hand, are limited due to th e relatively 
long exposure times requ i red fo r detections (e.g. lBeelen et alJ 
12006b IWang et al.1 l2008l 1201 Oh . Only about 10 (sub)mm 
quasars reported to date have detailed rest-frame FIR SED 
measurements. 

The Hersch el Space Observatory 

( Pilbr atTet al.1 1201 Oh has opened a new window 

dsPIRE: 250. 350. 500 /im. Griffin et all 1201 Oh to directly 
study the rest-frame FIR properties for quasars with moderate 
redshift s (z ~ 1.5). The Herschel Multi-tiered Extrag alactic 
Survey dHerMES. Oliver et alj|2010t lOliver et al.ll201 ll) cov- 
ers ~ 70 deg 2 with rich multi-wavelength data. Along with 
other Herschel surveys, rest-frame 30-300 /im emissions, 
likely due to cold dust in quasars and o t her AGNs, have been 
detected (jH atziminaoglo uet al.l 1201 Ot iLeipski et alJ 120101 : 
ISerjeantetal.ll2010h . 

In this paper, we report 32 Herschel SPIRE detections of 
24 /im flux-limited quasars in the HerMES Lockman Hole 
field. This SPIRE-detected sample constitute at least 10% 
of the 24 /im selected quasar sample in this field, and allow 
construction of the complete FIR SEDs for these broad-line 
quasars at z ~ 1.5, which triples the size of (sub)mm observed 
quasars that have detailed FIR SEDs. Throughout the pa- 
per, we assume a concordance cosmology with Hq=70 km/s 
Mpc -1 , ft M =0.3, and n A =0.7. 

2. MIPS 24/iM-SELECTED QUASARS AND THEIR FIR 
COUNTERPARTS 

The quasars used in this paper are from a 24 /im flux-limited 
sampl e in the Spitzer Wide-area In fraRed Extragalactic Sur- 
vey (ISWIRE: Lonsdale et al.ll2003l) . In the Lockman Hole- 
SWIRE (LHS) field, we selected targets that satisfy MIPS 
S 2 4 > 0.4 mJy (~ Sex), and 94% of the flux-limited sample 
also satisfy SDSS r AB < 22.5. In 2009, Huang et al (in prep.) 
performed a spectroscopi c survey of ~ 3 000 such 24 /im tar- 
gets with HECTOSPEC dFabricant et al.ll2005h on the Multi- 
ple Mirror Telescope (MMT), with an effective coverage of 
~ 8 deg 2 . 93% of these objects have re liable redshifts. SPSS 
objec ts with existing spectroscopic z (Hatziminaoglo uet alJ 
2008) that satisfy the same flux-limits were later added to 
the 24 /im flux-limited sample, which increased the spectro- 
scopic completeness to ~ 70%. Broad line quasars were 
then selected , where Mgll or CIV l ine width has a FWHM > 
1000km s" 1 (ISchneider et al.ll2007h . The final sample of 326 
24 /im-selected sources includes 210 MMT and 116 SDSS 
quasars. 

We matched these 326 quasa rs to the HerMES SPIRE cross- 
identification (XID) catalog (Ros eboom et alJ L2010I) . The 
XID catalog used SWIRE MIPS 24 /im positions to mini- 
mize the source blending effects due to large beam sizes (18" 
FWHM for SPIRE 250 /im images), and has a completeness 
of ~ 80% at S250 = 20mJy. Among the 326 24 /im-selected 
quasars, there are 41 SPIRE detections with S/N > 5, of which 
three were detected at 350 or 500 /im only. We dropped four 
sources whose SPIRE 250 /im beam covers two 24 /im coun- 
terparts (e.g. Fig. Q] bottom). We also excluded the five z 
< 0.5 objects because their rest-frame FIR data points do not 
constrain the SED fitting. The final sample consists of 32 
quasars (20 MMT and 12 SDSS objects) at 0.50 < z < 3.54, 
with a median z of 1.55 (Fig. [2] inset). This corresponds to a 
10% detection rate. Since 29 of them were SPIRE 250 /im de- 



tected, the 32 quasars used in this paper are hereafter referred 
to as FIR-detected quasars (Fig. Q] top and middle panels). 
Twenty-seven sources have at least one 350 or 500 /im detec- 
tion (> 3a), and 16 sources also have SWIRE MIPS 70 or 
160 /im detections. 

3 . SPECTRAL ENERGY DISTRIBUTION 

We constructed the rest-frame SEDs for the 32 FIR- 
detected quasars from the UV to the FIR bands. The LHS field 
was covered by the Galaxy Evolution Explorer (GALEX) in 
the ultra-violet (FUV, NUV), SDSS in the optical (u, g, r, 
i, and z), the UKIRT Infrared Deep Sky Survey (UKIDSS) 
in the near infrared (NIR; J,H,K), and the SWIRE survey in 
the mid-infrared (IRAC at 3.6, 4.5, 5.8, and 8.0 /im; MIPS 
at 24, 70, and 160 /im). A Chandra X-ray survey covered 
a small fraction of the LHS field (0.7 deg 2 ), and only one 
source (LHS-S 119) was d etected out of the t hree FIR-detected 
quasars within that area (Wilkes et al. 2009). Fig.Q]shows the 
stamp images for two FIR-detected quasars and one poorly 
matched quasar in the optical and infrared bands. 

The rest-frame FIR emissions of these quasars fall in 
the 30-300 um region, simi lar to those of ISO observed 
PG quasars (|Haas et al. 2003) and (sub)mm-detected quasars 
(e.g.lMcMahon et al.lll999b IWillott etal]|2003l; iPriddev et alJ 
I2003albb iRobson et al.l 12004b lOmont etaiT boo ll 120031 
I Carilli et al.l200lllBeelen et al.l2006tlWang et al l2007Ll2008l 
2010). Using the same method as was adopted in Beele n et alJ 
(2006) 19 , we derived the dust mass for the FIR-detected 
quasars. Md is in the order of 10 8 -10 9 M Q , similar to that 
of (sub)mm quasars. These values are 1~2 dex higher than 
values estimated for the PG quasars. Therefore, these FIR- 
detected quasars are associated with the (sub)mm detected 
quasars as a 'dust-rich quasar' population (Fig. |2}. 

In 3 1 out of the 32 quasars, the rest-frame SEDs exhibit FIR 
excess over the E94/R06 quasar templates (Sec. Q3 by 0.5 to 
2.3 dex at 90 /im, with an average of 1.4 dex. This suggests 
that the contribution from cool dust is present. We compared 
the mean SEDs for FIR-detected and the much larger sample 
of undetected quasars (Fig. 0). These mean SEDs were con- 
structed by combining individual rest-frame SEDs. We first 
converted the flux densities to luminosities for each object. 
After shifting their bandpasses to the rest frame, we normal- 
ized each SED at 1-5 /im to the R06 template. We then popu- 
lated a grid of points separated by 0.03 in log wavelength and 
linearly interpolated between the effective detections from the 
UV to the FIR. The mean luminosities with this gap-repaired 
photometry were then connected as the average SED (ma- 
genta curves). For FIR-undetected objects, their mean SED 
was compensated in the FIR with the stacked mean fluxes at 
250, 350 and 500 /im from the SPIRE images. These val- 
ues were estimated by first cutting out maps around the MIPS 
positions for individual sources (interpolation into sub-pixels 
was allowed for precisely centering), and then measured from 
the stacked map via a centered point-spread function (PSF) 
fitting. Errors associated with the mean fluxes were calculated 
using the bootstrap method. SEDs for both populations re- 
semble the R06/E94 templates in the optical and the NIR, and 
differ mainly in the FIR: the stacked S250 f° r FIR-undetected 
sources is ~ 8 mJy, about 4 times lower than the median S250 
for FIR-detected quasars (31.1 mJy). UV-optical reddening is 
common in both populations, being present in ~ 40% of the 

19 M d = S„ Dl 2 /(1 +z)ki(u)B(u, T d ), where k d (u) = k oju/upf is the dus t 
absorption coefficient. Here we used S250, and kj is from Alton et al. 12004). 
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FIR-d etected quasars. The reddening corrections are compli- 
cated (Hao et al. 2005) and beyond the scope of this paper, as 
we concentrate in the FIR. 

3.1. Modeling the FIR quasar SED 

For this work, we adopted a T-a-/3 model from Blain et al. 
(2003) to estimate the dust temperatures and quasar luminosi- 
ties. Different from a pure modified blackbody (MBB) model, 
which uses an exponential thermal function with emissivity 
index /3 to account for a single temperature dust component 20 ; 
in the T-a-f3 model, a power-law Wien tail (f u oc v~ a B v i) is 
introduced to the mid-IR SED to account for the warmer dust 
components (Fig. |4). This additional term is then matched 
to the MBB component at a transition point, where the two 
functions also have equal first order derivatives. This transi- 
tion w avelength varies fro m case to case. We adopted (3 = 2.0 
here (iPriddev et al.l2003al) 21 . At shorter wavelengths, we nor- 
malized the R06 template to each SED over the 1-5 /im range 
for reference. 

Out of the 32 quasars, four only have a single band rest- 
frame FIR detection. One (LHS-M020) of the four sources 
matches well to the R06 template in the optical and near in- 
frared bands. The quasar with the highest z, LHS-M268, also 
has an SED that is well described by the R06 template alone 
(within 0.2 dex in the FIR). So no fit was attempted for these 
five objects. 

We carried out T-a-j3 fits to the remaining 27 SEDs, of 
which 9 objects have all three SPIRE detections. Longward 
of the Lyman break (912 A), twenty-two quasars are well- 
defined by the R06 template and a T-a-/3 fit. All 12 SDSS 
quasars fall into this category. The remaining five less-well- 
defined SEDs show a strong stellar bump in the optical-NIR 
regime, presumably from the host galaxy, while their FIR 
SEDs are well described by a T-a-(3 fit. Fig. |4] presents the 
individual SED for each of the 32 FIR-detected quasars. The 
plots are labeled with their ID, redshift, fitted dust temperature 
and a, and follow an a decreasing order to match Table. Q] No 
host galaxy correction was applied. 

3.2. Dust temperatures and luminosities 

The fitted dust temperature (Tj) for these 27 quasars has a 
range of 18K-80K, with 87% of the sources in the range of 
25-60 K, and a median and mean of 29 K and 34 K (Fig. [6l 
right). The fitting error for Lj is less than 10% in ~ 70% 
of the cases. However, it is worth noting that the T& derived 
from T-a-P fit is on average 30% lower than that from a pure 
modified blackbody fit, where no power-law term is present. 
The FIR luminosities are similar within 3 a between these two 
fits. 

The fitted a for these 27 sources has a wide range of 
0.68 to 2.44 (Table. [TJ. Starbursts are found to have higher 
a values than normal sta r form ing galaxies and quasars, as 
was shown in iBlain et alj d2003l) . Different a values suggest 
different dust temperature compositions in individual quasar 
systems, and may be associated with different evolutionary 
stages: flatter slopes (a < 1.0) imply infrared SEDs with rel- 
atively stronger warmer dust emissions, likely heated directly 

20 f v oc v^B v ~[, where B„.x is the blackbody spectrum. 

21 A change in /3 from 2.0 to 1.5 do es no t significantly (< 3a) change 
the fitted and the Lfir (see also Sec. 13.21 . Since real dust may have a 
temperature distribution, the fitted temperature only applies to the cold dust 
component defined by the rest-frame FIR data, while the fitted a term indi- 
cates the relative strength of warm and hot dust. 



by the quasar; while steeper slopes a > 2.0 indicate a colder 
dust dominant infrared SED, similar to that of star forming 
galaxies (Fig. [6] left). Majority of the 27 quasars (~ 70%) has 
an a value of 1.0 < a < 2.0, possibly in a mixed condition 
between the two extremes. 

The FIR luminosities Lfir (40-300 /im) were estimated by 
integrating over the fitted S EDs, while the t otal infrared lu- 
minosities Lir (8-1000 /im) (lKennicuttll9 98) were integrated 
over the observed SEDs up to the redshifted MIPS 24 /im 
data, and over the fitted SEDs at longer wavelength. Lfir 
ranges from 10 to 10 135 L©, while Ltr has a range of 10 115 
to 10 14,3 L©, qualifying most of their host galaxies as ultra- or 
hyper-luminous infrared g alaxies (ULIRGs, Ltr > 10 12 L^; 
HyLIRGs, L m > 10 13 L©. ISaiTders & Mirabel Il996l) : in the 
27 quasars with a T-a-ft fit, there are 8 ULIRGs and 13 
HyLIRGs. Twenty-one of the 27 quasars have Lfir/Ltr > 0.2, 
and 16 show Lfir/Lir > 0.5, indicating major contribution 
from the FIR to the total Ltr. For FIR-undetected quasars, 
this ratio is ^0.3 based on their mean SED (Pig. [3]». We also 
estimated the 'big blue bump' luminosity, Lbbb (0. 1-0.4 /im), 
and the near-infrared luminosity, Lnir (2-10 /im), by integrat- 
ing over the observed SEDs for all quasars. The 2-10 jim 
range was chosen to minimize the stellar contribution, and 
to better represent the 'hot du st bump' emission likely to be 
directly heated by the quasar (Wang et al. 2008). These inte- 
gration ranges can be found in Fig. |3]as shaded regions. The 
parameters used and derived from the SEDs are summarized 
in Table. Q](See also Fig. [6]). Errors for Lj, a, Lfir and Lir are 
fitting errors only. The errors for Lbbb and Lnir contain only 
the photometric errors, and are not listed in Table. [T]because 
of their small values: mostly at the one percent level or less 
for Lbbb, and at most a few percent for Lnir. 

The FIR-detected and undetected quasars have similar Lbbb 
and Lnir values and distributions, with Lnir ^25% higher 
for FIR-detected quasars at higher redshifts (Fig. [7] left). 
The Lfir /Lnir ratio, however, shows an obvious excess over 
the E94 template range especially at Lnir < 10 13 L© (Fig.|7J 
right). We also plot the mean Lfir/Lnir ratios for FIR- 
undetected quasars in three Lnir bins. These values were de- 
rived from the mean SEDs at each Lnir bin in the 1 < z < 2 
range (covering the mean z of ~1.5), assuming ft = 2.0 and 
a = 1.6 (mean for FIR-detected quasars). These SEDs were 
completed in the FIR using the stacked fluxes from the SPIRE 
images for the relevant FIR-undetected sources. For com- 
parison, FIR-detected quasars in the same redshift range are 
marked with a cross in the center. 

It is clear that at the same redshift ( 1 < z < 2 in this case) and 
luminosity bins, FIR-detected quasars have higher Lfir/Lnir 
ratios than FIR-undetected objects, indicating differences in 
ALj, 7^, or AGN activities. This ratio seems to decrease at 
higher Lnir regardless of redshift, though the large scatter at 
Lnir < 10 13 Lq prevents determination of a tight correlation. 
The quasars' FIR properties cannot be predicted from shorter 
wavelengths, given the similar SED shapes and luminosities 
between FIR-detected and undetected quasars at optical and 
near-infrared regions (Fig. |7] left). Considering this similarity 
and the common UV extinctions, no bolometric luminosity 
was calculated for the FIR-detected quasars. 

4. DISCUSSION 

We constructed a sample of 32 quasars (0.5 < z < 3.6) with 
Herschel SPIRE detections, which triples the size of (sub)mm 
observed quasars that have detailed FIR SEDs. These FIR- 
detected quasars, as well as some FIR-undetected quasars in 
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our sample (Pig. [3j, show broad line features and strong cold 
dust emissions simultaneously. This is inconsistent with the 
evolutionary scenario that naked q uasars are only seen when 
the dust has been blown out (e.g. lHaas et alj|2003l) . The dust 
d etect ed by Herschel show temperatures from 18 K to 80 K 
( 33.2I >. with 87% of the sources in a range of 25-60 K, sim- 
ilar to local a nd z ~ 2 starburst g a laxies and UL IRGs (20 K 
< Th <60K) (ICalzettietalJl200aiMagdiretan [2010). The 
median and mean 7d in our sample are 29 K and 34 K, re- 
spectively. Since the FIR emissions for these FIR-detected 
and the (sub)mm detected quasars (e .g. lOmont etalJl200U 
l200l ICarUli et al.ll200H iBeelen et al.l 120061) both fall in the 
rest-frame 30-300 /im region, we associated these two popu- 
lations as a 'dust-rich quasar' population (Fig. |2). Estimated 
dust mass confir med this connectio n (See also Sec. |3). ISO- 
detected quasars (lHaas et al.ll2003l) are possibly at the low-z 
end of this population, though their could to be 1 —2 dex 
lower. A common assumption is that the dust-rich quasars 
are in transition between optically obscured and unobscured 
quasar phases, but near-infrared/(sub)mm spectroscopy and 
high resolution images are needed to test this. 

The FIR-detected and undetected quasars have similar 
mean SEDs (Fig. [3), redshift (Fig. [2] inset) and luminosity 
distributions (Fig. [7] left) at shorter wavelengths (rest-frame 
0.3-20 /im). Lnir is ~ 25% higher for the FIR-detected pop- 
ulation, likely affected by their FIR excess. The lack of 
correlation between properties at the FIR and shorter wave- 
lengths tests the widely used conversion factor of 7 between 
L(510 A) and Lboi for AGN and quasars (Shemm er et al.l 
120041) . For these FIR-detected broad line quasars (~ 10% 
of the flux-limited sample), if the rest-frame FIR emission is 
mainly due to quasar heating, this factor should be modified 
to 9—20 based on the Lfir/Ltr ratio. Whether this applies 
to the flux-limited quasar population in general remains to be 
investigated. 

At 1 < z < 2, the Lfir/Lnir ratios for FIR-detected 
quasars are on average 2x higher than FIR-undetected ones. 
Lfir/^nir and Lnir seem to be anti-correlated (Fig. [7] right). 
This trend is also observed for the overall population de- 
spite the large scatter at Lnir < 10 13 Lq: the relative strength 
of Lfir — com monly associated with star formation (e.g. 
iLutz et ai] l2007). decreases at higher Lnir, indicator of warm 
dust partially or mostly heated by AGN. This trend is consis- 
tent with the assumption that star formation is suppre ssed by 
the presence of a powerful AGN dHopkins et al. 2006). 

Both star forming galaxies and AGNs may contribute to the 
rest-frame FIR emission at various redshifts. If we attribute 
the FIR luminosity to star formation 22 , as was adopted in 



some previous studies (e.g. lEvans et : al.ll2.006l; iRiechers et al] 
l20MlNetzeret^l2lmlWang et al.ll2010Ll201 ll) . about 40% 
of the sample will require a star formation rate (SFR) > 
lOOOM^yr" 1 , si milar to that of submillimeter galaxies (e.g. 
ILutz et all 120081) . However, for the high luminosity end 



(Lnir >3x 10 13 L q , Fig. [7]), with only one source, the statis- 
tics are insufficient to prove whether there is a similar propor- 
tion of starburst dominated quas ars, as was found fo r (sub)mm 
observed quasars (i.e. -20-30%, Wan g et alJl2008h . 

On the other hand, for some quasars, the SFR derived from 
Lfir reaches 5000 M yr -1 , which is unlikely and probably 
implies that part of the FIR emission is powered by AGN. 
In Fig.|7](right), we found that for the more luminous quasars 
(Lnir > 5 x 10 12 L Q ), the Lfir/Lnir ratio falls in the normal 
quasar range, suggesting pure quasar heating. Resolved CO 
and PAH emission from NIR/mm spectroscopy and high reso- 
lution imaging showing dust distribution will provide a better 
estimate of the relative contributions from AGN and starbursts 
for these dust-rich quasars. 
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22 Eq(4) in lKennicutll fl993) : SFR (M yr 4 ) = 1.7 X 10' 10 Lir(L o ). Here 
we used ^Jr 300 instead of Lj^ 1000 for a more conservative estimate that re- 



duces further contamination from the AGN. 
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FIG. 1 . — Multi-wavelength images for two FIR-detected (top and middle) and one confused (bottom) quasars. The red circles in all three bands indicate the 
optical positions and have a 5 " radius. The sizes of the optical-mid-IR images are 50 " X 50 ", corresponding to the yellow squares in the 250 fim images. The 
cyan contours in the 250 fim image mark the 3, 6, and 9 X <r. The FIR emission corresponds to an isolated 24 source in all the 32 FIR-detected quasars listed 
in Table. [T] The bottom panel shows one of the four confusing cases excluded from the 41 sources with 5cr SPIRE signal. 
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FIG. 2. — Relation between the FIR luminosities and absolute B band magnitudes (left panel), and redshifts (right panel) for Herschel FIR-detected quasars 
(red), ISO-detected PG quasars (green), and (sub)mm-detected quasars — olive and yellow stars for quasars with single band (sub)mm photo metry , purple and 
blue squares for those with multi-band (sub)mm data (for references, see Sec. [3). Lfir was estimated with a T-a-/3 method for all sources (Sec. 13. if . The shaded 
area in red plots the SPIRE selection function at S250 = 20mJy, and the range results from the assumed 20% flux error (5 <r). This selection function was derived 
using fixed a = 1.6 and 7j = 34 K (mean values of the 27 objects with a T-a-/3 fit). Some sources are below the SPIRE selection region due to lower T$ or higher 
a. For objects with onl y one photometric point (850 Atm or 1 .2 mm data), we adopted for the fit a = 1 .6, and T& = 43K — mean values for (sub)mm quasars with 
multi-wavelength data i Beelen et al. 2006; Wang et al. 2008, 2010). The inset shows the redshift distribution for the FIR-undetected quasars in our sample, with 
the 32 FIR-detected quasars in red hatching and scaled by 4 X for comparison. 
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FIG. 3. — Rest-frame mean SEDs for FIR-detected and undetected quasars (arbitrary units). The colored dots show the detections (normalized to R06 template 
at 1-5 /an): grey, GALEX; purple, SDSS; cyan, UKIDSS; green, SWIRE-IRAC; orange, SWIRE-MIPS; red, HerMES-SPIRE. In magenta are the mean SEDs 
for FIR-detected and undetected quasars derived from these normalized detections (see Sec. [5). Dashed curves in dark and light blue mark the relative positions 
of the E94 and R06 templates (normalized at 1 (im). For FIR-undetected quasars, the stacked mean fluxes from the SPIRE images are plotted as red squares, 
shifted to the rest-frame assuming z = 1.5 (mean for FIR-undetected objects), and errors in the x-axis show the redshift span of 0.2 ~ 4.0. The regions used for 
L\,bh, ^NIR, and ^-FIR integration are shaded in yellow. 
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FIG. 4. — Individual rest-frame SEDs for the FIR-detected quasars in our sample, with energy \L\ in Lq. The SEDs were constructed with existing data from 
the UV to the mid-infrared, complemented in the FIR by the HerMES-SPIRE observations (see Fig.[3]for details). The red curves plot a T-a-fi fit with /3 = 2.0, 
with fitted Ti and a given in the legend, and blue curves are the R06 template normalized at 1-5 fim. The parameters are summarized in Table. [T] The objects 
are arranged in an a decreasing order, and followed by quasars with special SED shapes. No host galaxy correction was applied. 
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FIG. 5. — Fig.|4] Continued. 
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FIG. 6. — Histograms of the fitted power-law slope a (left) and the dust temperature Lj (right) for the FIR-detected quasars. 
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FIG. 7. — (left) Ratio between the near-infrared and the 'big blue bump' luminosities (definitions in Sec. l3.2l and Fig. [4] left). Hatched area marks the E94 
predicted range (90% confidence, E94 was used since R06 template stops at rest-frame 100 fim). ( right) The relation between FIR and NIR luminosities for FIR- 
detected quasars. Eleven FIR-detected sources have a higher ratio than the E94 prediction. Significant higher Lfir/Lnir ratios are observed at Lnir < 10 13 Lq , 
whereas more luminous quasars (Lnir > 1O 13 L0) show relatively weaker FIR dust emission. Large black crosses are the mean Lfir/Lnir ratios for FIR- 
undetected quasars at 1 < z < 2 in each Lnir bin. FIR-detected quasars of the same redshift range are marked with a small cross at the center. 
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Parameters derived from the fitted SED. 
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NOTE. — (1): Object ID. (2): Spectroscopic redshifts dete rmine d from Hectospec and SDSS spectra. (3) (4): Dust temperature 
and power-law index derived from the 1-a-f} fit (see also Sec. l3.lt . (5) (6) (7) (8): Luminosities estimated from the SED (see also 
Sec. l3.2t . (9): Star formation rate estimated using Z.p IR ^ m following the Kennicutt law (1998), see Sec.0 footnote. (10): Dust 
mass calculated using S950, see Sec. [3] footnote. 

*: not corrected for host galaxy contamination, f : SED well defined by R06 template, no T-a-fi fit needed. 



